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ABSTRACT 

In the second paper of this series, we present results from cosmological simulations on the demographics of 
flyby interactions to gauge their potential impact on galaxy evolution. In a previous paper, we demonstrated 
that flybys - an interaction where two independent halos inter-penetrate but detach at a later time and do not 
merge - occur much more frequently than previously believed. In particular, we found that the frequency of 
flybys increases at low redshift and is comparable to or even greater than the frequency of mergers for halos 
> 1O'^/!“^M0. In this paper, we classify flybys according to their orbits and the level of perturbation exacted 
on both the halos involved. We And that the majority of flybys penetrate deeper than ^ Rhaif of the primary and 
have an initial relative speed ~ 1.6 x Vyir of the primary. The typical flyby mass-ratio is ^0.1 at high z for 
all halos, while at low z, massive primary halos undergo flybys with small secondary halos. We estimate the 
perturbation from the flyby on both the primary and the secondary and And that a typical flyby is mostly non- 
perturbative for the primary halo. However, since a massive primary experiences so many flybys at any given 
time, they are nearly continually a victim of a perturbative event. In particular, we And flybys that cause ~ 1% 
change in the binding energy of a primary halo occurs > 1 Gyr“' for halos > 1O^°/!“'M0 for z 4. Secondary 
halos, on the other hand, are highly perturbed by the typical encounter, experiencing a change in binding 
energy of nearly order unity. Our results imply that flybys can drive a significant part of galaxy transformation 
at moderate to lower redshifts (z < 4). We touch on implications for observational surveys, mass-to-light ratios, 
and galaxy assembly bias. 

Subject headings: cosmology: theory — cosmology: dark matter — cosmology: large-scale structure of uni¬ 
verse — galaxies: evolution — galaxies: halos — galaxies: interactions — methods: nu¬ 
merical 


1. INTRODUCTION 

Galaxy mergers drive galaxy evolution; they are a key 
mechanism by which galaxies grow and transform. The 
merger process is usually accompanied by a strong gas in¬ 
flow (??) triggering a central starburst and perhaps an ac¬ 
tive galactic nucleus (hereafter, AGN). Consequently, galaxy 
mergers have been studied quite extensively both numeri¬ 
cally (e.g., ??????) and observationally (e.g., ?????). Galaxy 
mergers are so successful in driving galaxy evolution simply 
because they strongly perturb the potential. Of course, weaker 
perturbations can cause change, too; it has long been thought 
that even an orbiting low mass satellite like the Large Magel¬ 
lanic Cloud can distort the underlying smooth galaxy poten¬ 
tial enough to warp the Milky Way (MW) disk (?). However, 
one entire class of galaxy interactions capable of causing such 
perturbations - flybys - has been largely ignored. 

Unlike galaxy mergers where two galaxies combine into 
one remnant, flybys occur when two independent galaxy halos 
inter-penetrate but detach at a later time; this can generate a 
rapid and large perturbation in each galaxy. We define a flyby 
as a fast interaction between two independent galaxies each 
within their own separate dark matter halos. Although some 
aspects of a galaxy flyby may be similar to the first wide pass 
of a galaxy merger, or to galaxy harassment in a cluster, the 
rapid and transient nature of a flyby places it in a dynamical 
class of its own. 

One of the first questions to ask is how common are galaxy 
flybys?. Curiously, though, this is a difficult question to an- 
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swer simply because previous theoretical work has been set 
up to search for mergers; indeed both the extended Press- 
Schechter framework hierarchically assembles a dark matter 
halo, and the merger trees gleaned from cosmological Nbody 
simulations are designed to track halo growth in the same way 
- solely through mergers and smooth accretion. 

Previous numerical works have shown the existence of fly¬ 
bys, or ‘backsplash’ galaxies, in galaxy clusters (e.g., ??) and 
isolated galaxies (??). Simulations predict that up to ~ 60% 
of the galaxies found between 1 -2Rvir could be ‘backsplash’ 
galaxies (???). Naturally, these galaxies are distinct from 
those that are infalling for the first time. Backsplash galax¬ 
ies are subject to ‘pre-processing’ and usually exhibit HI defi¬ 
ciency compared to infalling satellites. Candidate backsplash 
galaxies have been identified in galaxy clusters (???), as well 
as the held galaxies in the Local Group (?). Another theoreti¬ 
cal approach at finding flybys focused on identifying subhalo 
orbits that extend far outside the primary virial radius (???). 
? looked at a suite of five isolated Milky-Way sized dark mat¬ 
ter halos in a cosmological simulation and found that roughly 
half of the subhalos that were once associated with the host 
are now located outside the virial radius. Broadly speaking, 
flyby halos constitute ^ 50% of all halos on scales of 1 -3 Rvh . 

Recently, we developed and tested a method to identify both 
mergers and flybys between dark matter halos in cosmological 
simulations (?). We constructed a full ‘halo interaction net¬ 
work’ that assesses the past dynamical history of any given 
halo. With this new too|^ we made the first census of fly¬ 
bys, and we found that they are surprisingly common. In fact 

^ available publicly at https://bitbucket.org/manodeep/ 
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at z < 2 , flybys occur more often than mergers for all halos 

> 1O"/!-'M0. 

We reiterate that majority of flybys are one-time events, not 
merely the first orbit of a merger - here, the disks (if present) 
will not physically overlap at any point; they are also funda¬ 
mentally different than galaxy harassment. More cosmolog¬ 
ical simulations are underway to build up a better statistical 
sample and to expand the dynamic range in halo mass, but 
the emerging picture is that current studies focusing on solely 
merging galaxies may completely miss this common type of 
galaxy interaction. 

Now that we have evidence that flybys are rampant, we 
need to know if neglecting these flybys matters. While dark 
matter halo flybys do not significantly affect the halo mass 
function, these fast, transient events may make their mark 
on galaxy structure and kinematics. Over the years, flybys 
have been invoked in a heuristic way to explain various mor¬ 
phological features and transformations, such as the evolution 
from spiral to SO galaxies in group environments (?), or the 
excitation of spiral arms in the galactic disk (e.g., ?). Re¬ 
cently, flybys have been proposed as a way to provide the ex¬ 
tra starbursts needed to explain the steep abundance gradient 
found in massive ellipticals (?), because it is widely-held that 
these galaxies are too a-enhanced to be explained by merg¬ 
ers or galaxy harassment alone (??). From the purview of 
linear perturbation theory, a flyby imparts an impulse to the 
galaxy (?), and in principle this can exact a morphological 
and kinematic change in both the ‘victim’ and the ‘perturber’ 
halo. In fact, high resolution isolated simulations have borne 
this out, showing that flybys can induce a strong and long- 
lived bar (?), can spin up the dark matter halo (?), and can 
generate an S-shaped warp as well (?). 

In this paper, we aim to quantify the flyby perturbation and 
the frequency of disruptive flybys as a function of halo mass 
and redshift. In Section]^ we briefly describe the simulation 
details and the interaction network (see ? for further details), 
in Section |3] we show the results and Sections |4] and |5] we 
discuss the implications of our results and future steps. In 
addition, we outline the details of the analytic approximation 
used to estimate the perturbation on the halo in Appendix [A| 

2 . METHODS 

We use a high-resolution dark-matter simulation with 
1024^ particles in a 50 /r“'Mpc box with Ob = 0.044, 0dm = 
0.214, Oa = 0.742, ag = 0.796,/i = 0.719, n, = 0.96, consis¬ 
tent with WMAP-5 cosmology (?), where the symbols have 
their usual meaning. The initial particle distribution is ob¬ 
tained from a Zel’dovich linear approximation at a starting 
redshift of z = 249 and is evolved using the adaptive tree- 
code, GADGET-2 (??). The dark matter particles have a fixed 
gravitational softening length of 2.4 co-moving /!“*kpc at all 
epochs. We store 105 snapshots spaced logarithmically in 
scale-factor, a = 1/(1 -I-z), from z = 20 to 0. This translates 
into a timing resolution < 50 Myr for z > 3 and ~ 150 Myr for 
z < 3; in SHI2, we showed that a high snapshot cadence was 
essential in fully detecting the flybys in the volume. Since the 
fundamental mode goes non-linear at z = 0, we only present 
results up to z = 1 where the 50 /j“'Mpc box is still a represen¬ 
tative cosmological volume. 

To identify halos, we use a canonical linking length b = 0.2 
(~ 10 /!“'kpc) to find the Friend-of-Friends (EOF) halos (?); 
subhalos, with a minimum of 20 particles (~ 1 O*^/!“'M 0 ), are 
identified using the SUBFIND algorithm (?). The SUBFIND al¬ 


gorithm identifies subhalos as bound structures around a den¬ 
sity maxima; the background FOF obtained at the end is com¬ 
prised of particles that are bound in the general potential of 
the FOF, but not bound to any subhalo. For the remainder of 
this paper, all references to the FOF halo will mean this bound 
background halo. For details on the SUBFIND algorithm, we 
refer the reader to ?. We used an overdensity, pA = 200pcrit(2) 
to compute the virial radius of halos: 


Rv 


/3M™ 

V^ttpa 


1/3 


( 1 ) 


Here, we used the bound mass of the halo reported by 
SUBFIND as Mvir. Another commonly used overdensity is, 
PA.vir ~ (187r^-l-823c-393:^) x Pcrit(^)(?), where x= 0„,(z)-1. 
For ACDM cosmology at z = 0, pA.vir ~ 96 pcUt- For compar¬ 
ison, the virial radii obtained at z ~ 0 by using pA.vir would 
be roughly (200/96)'/^ ~ 1.27 times larger than our adopted 
definition. At higher redshifts, pcUt ~ Pmatter, and there is vir¬ 
tually no difference between our adopted pA and pA.vir- 


2.1. Characterizing flyby orbits 

Our next step is to identify all interactions between halos 
of any type in the simulation. This done with our halo in¬ 
teraction network algorithm described in SHI2. For every 
flyby we record the initial separation, rsep, a minimum sep¬ 
aration (impact parameter), b, between the halo centers dur¬ 
ing the encounter and the initial relative velocity between the 
halo centers, Vsep, at the beginning of the encounter. These 
three quantities, along with the two halo masses and the ini¬ 
tial redshift are sufficient to crudely specify the interaction. 
In other words, idealized galaxy flyby simulations can be run 
with these 6 parameters specified. We use physical values for 
all the quantities to facilitate such idealized setups. We track 
the orbit of the secondary halo during the flyby and mark the 
minimum separation between the primary and secondary cen¬ 
ters as b. We also note the primary half-mass radius (deter¬ 
mined from the particle distribution) at the time of this mini¬ 
mum separation. However, since b is the minimum separation 
found during snapshot outputs, it is really an upper limit on 
the true minimum separation between the two halos. To better 
estimate the true impact parameter, we numerically integrated 
the orbit of two analytic dark matter halos with the same ini¬ 
tial total mass, separation, and velocity as our target halos. 
To simplify, we assumed the halo concentrations were given 
by (?). As we integrate the orbital trajectory, we include the 
effect of tidal stripping to the Roche lobe, as well as Chan¬ 
drasekhar dynamical friction with a Coulomb logarithm that 
goes as log(l -l-(Mprim/44sec)^)- In practice, the mean impact 
parameter from the simulations was only 10% times larger 
than from integrating the orbit, but there was a non-negligible 
fraction of impact parameters that were larger by order unity 
- these encounters were typically highly radial and the likeli¬ 
hood of the simulation snapshot catching it at its true pericen- 
ter is small. We will use the b determined from this extrapo¬ 
lated trajectory for the rest of the paper. 


2.2. Perturbations from Flybys 

As a first estimate of the potential damage done by flybys, 
we follow the prescription set out by ?. Their work is based 
on linear perturbation theory to determine the change in to¬ 
tal and potential energy for spherical systems. The perturba¬ 
tion equations have different forms depending on whether the 
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halo experiences an internal or external encounter (see Ap¬ 
pendix!^. The generic flyby consists of a smaller secondary 
halo fallmg into a more massive primary halo and then detach¬ 
ing at some later time. In the reference frame of the primary, 
it is a victim of an internal flyby; for the secondary halo, how¬ 
ever, the center of the primary halo always remains outside 
Rvir so it operates as an external flyby. Therefore, we use the 
external and internal flyby perturbation from ? to calculate the 
perturbations on the secondary and primary halo respectively. 
We combined the internal and external flyby perturbations fits 
into one equation: 


AEi 


M2 

O.lMi 


exp(-0.5/3) exp(/3fe')V: 


sep 


( 2 ) 


where we replaced K with exp(-0.5/3) x expiPb'), and we 
used a = 1.95, B_= -1.0. This provides a continuous fit for 
AE jE in Eqn. A4 while preserving the energy changes from 


the two equations. Essentially, we used an exponential cut¬ 
off from the internal flyby perturbation calculation and ex¬ 
tended the prescription to b > 2 Rhaif.primary Since this is 
somewhat of an ad-hoc implementation, the exact perturba¬ 
tion will not be accurate; however, in this paper we are only 


concerned with relative estimates and these equations should 
be sufficient for our purposes. In contrast, ? performed high- 
resolution simulations of galaxy flybys and directly measured 
the perturbed modes and their magnitudes; they found that 
the 0.1 mass ratio secondary during a deep flyby could excite 
a perturbation in the m = 2 mode as high as 20% (with larger 
perturbation amplitudes for lower-order modes), resulting in 
a distinct and long-lived bar. 

3. RESULTS 

In ?, we showed that flybys occur frequently for massive 
halos (> 10 "/!-'Mo) for z < 2. Here, we show that flybys 
are potentially an important player in galaxy evolution. 

Eirst, we compare the frequency of flybys vs mergers from 
a more observational perspective. Table lists the fraction 
of flybys and mergers for various host halos at a range of 
redshifts and relative velocities. It is striking that for z < 2, 
the potential flyby contamination in a galaxy-pair survey is 
as much as ~ 85%, and at minimum, the contamination is 
roughly 1/3. It appears that galaxy surveys designed to probe 
the effects of mergers on star formation, AGN activity, or mor¬ 
phology could be picking up a significant fraction of flybys as 
well, particularly for Milky Way-mass halos. In this case, fly¬ 
bys do ‘matter’ because they are an abundant contaminant. 


TABLE 1 This table lists the number of flybys and mergers for various redshift, primary 
mass, and different velocity offsets for galaxy pairs. Column 1 shows the redshift range 
for the beginning of the interaction, column 2 shows the virial velocity of the primary 
halo, column 3 shows the relative physical velocity of the secondary halo (i.e., a velocity 
offset), column 4 and column 5 show the number of flybys and mergers respectively while 
column 6 shows the percentage of flybys ( = Column4/(Column4 + Column5)). The 
mass-ranges corresponding to Column 2 are 10* * — 10*^ and 

10*^ —4 X 10 '^ . In principle, this table shows the fraction of interactions that will 

be flybys from an initial halo pair with a given primary virial velocity and a velocity 
offset. For z = 1—2, the flyby rate is at least comparable to mergers for all halo masses. 


redshift range 

Primary Vvir 

V infall 

nFlybys 

nMergers 

Flyby fraction 

H 

[km/s] 

[km/s] 

H 

[-] 

[%] 

1.0-2.0 

85-190 

170-380 

255-570 

340-760 

1789 

200 

32 

1082 

66 

5 

62 

75 

86 


190-645 

380-1290 

570-1935 

911 

29 

1256 

60 

42 

33 

2.0-3.0 

100-220 

220-750 

200^40 

300-660 

440-1500 

660-2250 

692 

37 

176 

9 

1022 

68 

567 

72 

40 

35 

24 

11 

3.0-4.0 

115-250 

250-845 

230-500 

345-750 

500-1690 

750-2535 

155 

3 

7 

0 

541 

19 

112 

15 

22 

14 

6 

0 


3.1. Two classes of flybys 

A flyby will be perturbative when the halos penetrate deeply 
before detaching; in addition, if the flyby is relatively slow, 
then the secondary halo particles spend more time near the 
resonances - thereby perturbing the primary even more. Thus, 
highly perturbative flybys would have to be highly penetrat¬ 
ing and relatively slow. In the next sections, we will dissect 
the orbital characteristics of the flybys in the simulation and 


examine their relative importance over a range of halo masses 
and redshift. 

Eirst, we identify the infall characteristics of flybys. Eor a 
flyby to occur, the secondary halo has to overcome the pri¬ 
mary halo’s gravitational potential; therefore, the initial rel¬ 
ative velocity must be high compared to the primary halo 
escape velocity. A priori, therefore, we expect that flybys 
will have an higher initial relative velocity compared to merg¬ 
ers. Nominally, an unbound interaction should have Vsep > 
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Fig. 1.— Probability density of infall velocities (scaled by the pri¬ 
mary V»ir) for mergers (blue, dash-dotted line) and flybys (red, solid 
line). The merger infall velocity peaks at ~ 1.2Vvir while the flyby 
infall velocity peaks at a higher value of ~ l.bVvir. Formally, the 
escape speed for the interaction is v/2Vvir. Hence, most of the fly¬ 
bys have infall velocities much larger than the escape speed of the 
primary halo. A two-sided K-S test shows that merger and flyby 
infall velocity distribution come from different populations (p-value 
for the null hypothesis that they are the same distribution is < 10^*). 

Vesc,primary = V^Vviri however, mass-loss and dynamical fric¬ 
tion effectively increase this minimum infall velocity. Previ¬ 
ous work looking at the distribution of subhalo infall veloc¬ 
ities has found that typical Vinfaii Vyir and only a few per¬ 
cent of the orbits are truly unbound (e.g. ??) In Fig. [T] we 
plot the histogram of the relative velocity scaled by the pri¬ 
mary Vvir for both mergers and flybys. We note that flybys 
and mergers are indistinguishable in a plot for the ratio of rsep 
to Rvii - Since we classify flybys based on past/future behav¬ 
ior depending on the velocities of the halos, there is no dis¬ 
tinction between mergers and ‘grazing’ flybys in real space. 
As expected, the infall velocity of mergers peaks at a lower 
(bound) velocity, ~ 1.2Vvir compared to the peak flyby infall 
velocity of ~ l.hVvir. We And that roughly 1% and 10% of 
the flybys occur with Vsep < v/2 and 2 times the primary Vyir 
respectively. For secondary halos, about 50% of the flybys 
occur at relative velocities > 5 x Vyir, while only ^ 30% of 
the mergers occur at Vsep ^5 x Vyir. 

Fig.[2demonstrates that flybys have a higher infall velocity 
compared to mergers but the initial separation between cen¬ 
ters for both mergers and flybys is approximately the sum of 
the virial radii of the two halos. So, the only way to obser- 
vationally distinguish between a flyby and a merger would be 
through the 3-d relative velocity, which makes it problematic 
to disentangle mergers from flybys in a survey. 

We can also characterize the distribution of the flyby im¬ 
pact parameter. In Fig. we show the probability density of 
b, scaled by primary Rhaif at the time of minimum separation, 
for three primary halo mass ranges and three mass-ratio cuts. 
We calculate the probability density with an adaptive kernel 
density estimator (?) with a Epanechnikov kernel. Combining 
the impact parameter results with the flyby velocity distribu¬ 
tion above, it is clear that there are two classes of flybys: one 
that rapidly delves into the primary core, and one that slowly 
skirts the outer primary halo. For massive halos, the most 
common flyby is a deep one, with few skirting encounters. In 
general, major flybys are of the slow and shallow type, while 



^ / ^half,pri 


Fig. 2.— The probability density of flyby impact parameters for 
different mass-ratios and primary masses. The top left panel shows 
all primary masses > IO^A^'Mq, the bottom left shows 10^- 
10*’ while the bottom right panel is for all primary masses 

above > 10*'The (red) solid lines shows all flybys with 
q> lO^”*, the (blue) dash-dotted line shows ^ = 10“'* -10“’, while the 
(golden yellow) dashed line shows major flybys with q = 10“’ - lO”. 
The two arrows show the 25% and 75% quartile for all the flybys 
(i.e., q > 10^) in that panel. From the plots, we can see that the typ¬ 
ical b for a flyby is ~ primary Rhaif for mass-ratio q< 10“’. The fly¬ 
bys with ^>0.1, almost inevitably are distant. We postulate that this 
is an outcome of efficient dynamical friction for close-to equal-mass 
interactions, where the subhalo loses its kinetic energy and cannot 
escape the main halo potential for deep orbits. Therefore, a typical 
‘major’ flybys occurs with b > 3 Rhaif and is essentially a grazing 
encounter. 


minor flybys penetrate deeper into the primary halo potential. 

3.2. How long do typical flybys last? 

We have so far looked at the distribution of b and Vinfaii for 
flybys. Now, another important orbital characteristic of fly¬ 
bys will be the typical duration. We remind the reader that a 
flyby is characterized by a ‘main halo—J^subhalo—J^main halo’ 
transition, and dynamically this means that the secondary halo 
must fall in to the primary halo and climb out of the potential 
well. This implies that the typical flyby duration is 0{tcross)- 
In fact, in ?, we used a minimum duration of half a crossing 
time to distinguish ‘true’ flybys from artificial halo stitching 
due to the FOF algorithm. We choose the crossing time at the 
beginning of the flyby, tcmss, to express the flyby durations. 
By selection, all flybys must last at least 0.5 fcrou- We note 
that tcross is independent of halo mass and evolves as 
with z. 

In Fig. we show the median flyby duration as a function 
of halo mass and redshift for both the primary and secondary 
halos. The duration is scaled by tcmss - where tcnss is eval¬ 
uated at the beginning of the flyby. From the figure, we can 
see that the median duration of a flyby is typically larger than 
two tcross- There is an redshift trend for the duration as well 
- high z flybys are shorter compared to the low z ones. For 
z < 5, both primary and secondary halos typically undergo 
flybys that last 2-3 tcross- However, since we can only detect 
the beginning and end of flybys at a snapshot (which goes as 
loga), we caution that the true flyby duration may be shorter. 
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Fig. 3.— Median flyby duration, in units of tcmss as a function of 
halo mass and redshift. Primary halos are plotted on the right and 
secondary halos are plotted on the left. The redshift and halo mass 
used for the plot is at the beginning of the flyby. Yellow/red means 
that for flyhys in that mass-z bin, the median flyhy will last >2-3 
tcross for that redshift. At z > 8 , flyhys are very fast and typically last 
one tcross (~ 50-70 Myrs) whereas at lower redshift, z < 5, flybys 
typically last 2-3 tcross (~ 1 - 2 Gyr). 

3.3. Distribution of mass ratios 

In any cosmological simulation, the halo number density 
is highest for the lowest resolvable halo. Consequently, the 
frequency of interactions increase with decreasing mass ratio, 
roughly as In Fig. we show the distribution of me¬ 

dian flyby mass ratios for primary and secondary halos. For 
primary halos, the trend is roughly independent of redshift, in 
other words, the mass ratio smoothly decreases with increas¬ 
ing halo mass. The lowest mass primary halos have flybys 
that are approximately equal-mass, while the twical flyby for 
the most massive primary halos is a o' ~ 10 ““*. n 

Note that interpreting this figure is not strai^tforward be¬ 
cause it also encodes the halo mass function. For example, 
for the massive secondary halos, e.g, lO" /!“'M 0 at z ~ 5, the 
only feasible primary halos are also of similar mass. Hence, 
the mass ratio for such flybys is close to 1. For a low-mass 
secondary halo, the typical interaction is q ^ 10 “* - the con¬ 
verse of its high-mass primary halo counterpart. Here we see 
that flybys very rarely occur among massive, equal-mass ha¬ 
los. While one reason is that any major interaction is always 
less likely compared to a minor one, we do still see major 
mergers in the simulation. So, naturally the question arises - 
why are there fewer major flybys? The answer comes from 
the velocity function of dark matter subhalos. In ACDM cos¬ 
mology, the velocity of infalling satellites drops of steeply, as 
exp(-(vjaf/vAojf)“), with a ~ 3 (?). In addition, the efficiency 
of dynamical friction and tidal stripping ensures that massive 
subhalos get stripped, lose significant amount of energy and 
angular momentum, and have a more difficult time making it 
out of the primary halo. This effect was pointed out by ?, who 
noticed the mass-bias in unorthodox orbits - where smaller 
subhalos can preferentially lie as far as 5 Rvii from the primary 
host halo. The initial velocity distribution and the preferential 
energy loss for massive halos combine to make major flybys 
comparatively rare. 

3.4. Mass loss from flybys 

The secondary halos are less massive than the primary halos by con¬ 
struction - therefore the mass ratio for the secondary halos are always greater 
than 1. 


Fig. 4.— Distribution of the median mass ratio of flybys as a func¬ 
tion of halo mass and redshift. Primary halos are plotted on the right 
while secondary halos are plotted on the left. Our formal halo de¬ 
tection limit is ~ lO*/i“*M 0 but we only use halos with at least 100 
particles, i.e., ~ 5 x halos. Thus, for the low mass primary 

halos, e.g., ~ lO*’/i“*M 0 , all interactions are resolution limited and 
appear roughly equal mass for all redshifts. In the same spirit, the 
typical interaction for a massive halo is with a ~ lO®/z“*M 0 and the 
smallest mass ratio can be as low as ~ 10“*. For a massive secondary 
halo, there are very few more massive halos and consequently, the 
typical interaction (if there is one), is more equal-mass. 

Now that we have discussed the demographics of flyby en¬ 
counters, we can turn our attention to the effect that flybys ex¬ 
act on the participating halos. This can address the question of 
whether flybys matter in transforming the halo structure and 
kinematics. One key effect is that the interaction can cause 
significant mass loss in the intruder halo. This is quantified in 
Fig. 0 which plots the fractional mass loss of the secondary 
halo as a function of the normalized impact parameter. This 
draws out two essential points: first, deep encounters strip the 
secondary halo mass by of order 50%, and second, the graz¬ 
ing encounters may sometimes gain mass as they travel along 
the equipotential contour of the primary halo. From this re¬ 
sult, we can already answer our original question of do flybys 
matter? with a definitive yes, because they dramatically alter 
the mass (and therefore the mass-to-light ratio) of the intruder 
halo. 

3.5. Flyby-driven perturbations 

Taking the concept of flyby-driven halo transformation a 
step further, we quantify the frequency of perturbative flybys. 
Broadly speaking, the perturbation from a flyby depends on 
the mass ratio, the impact parameter and the relative velocity. 
The usual intuitive understanding of interactions apply - more 
massive secondaries cause more damage, as do smaller impact 
parameters and lower relative velocities. 

With all the essential pieces now in hand, we will look at 
the perturbations from flybys - these perturbations are esti¬ 
mated using Eqn.|^ Looking at the functional form of Eqn.l^ 
the most important bit is the impact parameter, b which de^ 
termines if AE/E is characterized as an internal or exter¬ 
nal flyby. In Eig. we show the mean perturbation in 2-d 
cells of halo mass and redshift. Eor the secondary halos (left 
panel), we see that the mean AE jE during a flyby is of or¬ 
der unity and should strongly affect the internal structure of 
the secondary halc|^ However, for the primary halos (right 

^ tsEjE of order unity implies the assumed linear regime for deriving 
Eqn.j^is no longer valid 
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Fig. 5.— Mass loss of the secondary (1„/„//) as a function 
of the impact parameter (scaled by the host Rvir). There are two 
distinct peaks in the distribution: i) deep flybys and ii) grazing flybys. 
Grazing flybys typically do not lose much of their halo mass while 
the deep flybys can lose ~ 50% of the infall halo mass. 


panel), the mean AE/E from a flyby is essentially neglige- 
ble. Thus, for massive primary halos, a typical flyby is not 
going to change the overall structure at all. 

However, since the big question we want to know is if^y- 
bys matter, we need to know if any, and not just the typical 
flyby strongly affects a halo/galaxy. In Fig. we show the 
event rate of perturbative flybys, defined as AE/E > 10“^, 
as a function of halo mass and redshift. On the left panel, 
we show the event rate for ‘big’ flybys for secondary ha¬ 
los. The average event rate for secondaries increases with de¬ 
creasing redshift; with typical rates ~ 0.01 per Gyr at z > 4, 
and ~ 0.3 per Gyr for 1 ^ z ^ 4. The event rate resembles 
the distribution of durations (see Fig. and shows almost 
no trend with halo mass. The similarity with the duration is 
likely from the fact that deeper flybys take longer to com¬ 
plete, and these deeper flybys produce larger perturbations. 
For the primary halos, there is a strong trend with halo mass 
- with the event rate increasing by almost an order of mag¬ 
nitude between 1 O'°/!“'M 0 and lO"/i“'M 0 halos. All halos 
above 1 O"/!“'M 0 have at least a few ‘big’ flybys every Gyr. 
Thus, both primary and secondary halos have perturbative fly¬ 
bys frequently enough that the halo internal structure will be 
modified (??). 

In Fig. we saw that the typical flyby was highly pertur¬ 
bative for all secondary halos but in Fig. we see that the 
typical secondary halo only experiences ^1 highly perturba¬ 
tive flyby per Gyr. Similarly, in Fig. we saw that the typi¬ 
cal flybys don’t significantly perturb the primary, whereas in 
Fig. l7j we see that primary halos undergo a large number (few 
to l(r s) of highly perturbative flybys per Gyr. This apparent 
discrepancy between Fig. and Fig. l7j can be explained in 
terms of the hierarchical nature of halo interactions, and the 
perturbation calculation (Eqn.[^ itself. Secondaries don’t al¬ 
ways have flybys, but when they do, it is a highly perturbative 
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Fig. 6.— The mean of AE/E as a function of halo mass and red¬ 
shift. The primary mass is plotted on the right and the secondary 
mass on the left. We show the mean (and not the median) here to 
have a more representative measure of the perturbation on the ha¬ 
los. It is clear that the mean perturbation from flybys is much more 
damaging for the intruder halos than the primary halos. 
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Fig. 7.— The rate of big flybys (AE/E > 10“^) as a function of 
halo mass and redshift. The rate is plotted on the right for primary 
halos and on the left for secondary halos. The color shows the log¬ 
arithm of the flyby rate per halo per Gyr. Even though in Fig. we 
saw that the mean AE/E is small for massive primary halos, these 
massive halos undergo so many flybys that there is always some big 
perturbative flyby that occurs at any given z < 3. For instance, all 
halos > lO" at z < 3, have a perturbative flyby rate > 1 per 

Gyr. Such frequent disturbances will influence the evolution of all 
such massive halos. 



one. Averaged over the entire population of secondary ha¬ 
los, the event rate of highly perturbative flybys is small. On 
the other hand, massive primary halos undergo a lot of fly¬ 
bys; however, most of these are with small secondaries and 
the term in Eqn. reduces the typical AE/E drastically. 
Thus, the typical interaction strength for the primary halo is 
actually very small. But, in this large number of tiny flybys, 
there are a few interactions where the secondary gets close 
to the center of the primary halo. These flybys, with small 
impact parameters, produce really large AE/E and drive the 
high event rate for massive primaries. 

4. DISCUSSION 

Elybys and mergers have similar separation distributions, 
differing only by the relative velocities between the halos. 
One obvious implication of this similarity in separation is that 
it is very difficult to distinguish flybys from mergers in a pho¬ 
tometric survey, and even adding radial velocities will only 
help in a statistical sense. We should expect that merger rates 
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derived from pair counts will be contaminated with flybys. 
However, since flybys typically occur with low mass-ratios, 
this contamination affects minor merger rates much more than 
major merger rates. 

One consistent picture that emerges from flybys is that the 
secondary halos typically lose ^ 50% of their infall mass, and 
are highly perturbed as well. Previous work has shown that 
this level of perturbation drives bars and warps, and spins up 
the halo. Another potential side-effect of a flyby could be the 
shut-off in smooth gas accretion and consequent quenching of 
star formation in the secondary as it delves into the primary 
halo. 

In addition to these potential changes to internal galactic 
structure, flybys also have an effect on the overall galaxy dis¬ 
tribution. Flyby-induced mass loss, while keeping the bary- 
onic component of the secondary galaxy relatively intact, sys¬ 
tematically decreases its mass-to-light ratio (T*) at fixed lu¬ 
minosity. Observationally, these would imply a higher mass- 
to-light ratio (T*) and redder colors for a fraction of galaxies 
that are near other more massive galaxies. Since flybys are hi¬ 
erarchical, i.e., more massive halos have more flybys (?), the 
effect of flybys would be more prominent in galaxy groups 
and clusters. Such an excess in the red-fraction has been re¬ 
ported (?), and more recently ? showed that that the excess 
in red-fraction could be explained by assuming flyby galaxies 
evolved like satellites even outside the host. 

Another effect of tidal stripping during a flyby shows up as 
halo assembly bias. Halo assembly bias means that the clus¬ 
tering of halos depends on a hidden parameter beyond the halo 
mass. Over the past decade, halo assembly bias has been con¬ 
sistently detected in simulations (e.g., ????) and some have 
invoked ‘ejected’ subhalos as an explanation. However, a far 
more important consequence of flybys would be galaxy as¬ 
sembly bias - where galaxy properties depend not just on the 
host halo mass but another secondary characteristic, such as 
environment. If flybys systematically reduce halo masses, a 
galaxy-mapping prescription based only on halo mass would 
erroneously place a lower mass, bluer galaxy within. Or, 
worse still, the halo may not even receive a galaxy. Com¬ 
paring this mock galaxy catalog with observations, the mock 
would show a dearth of red galaxies around massive systems, 
as well as a suppression of galaxy clustering. As emphasized 
recently (?), galaxy assembly bias can cause Halo Occupa¬ 
tion Distribution (HOD) models to converge on a statistically 
different parameter space. Galaxy assembly bias from flybys 
would require an overhaul of the models that place galaxies 
into halos and statistically model galaxy formation physics, 
such as the HOD (?). Similar concerns apply for the Condi¬ 
tional Luminosity Function (CLF) (e.g., ?), as well. 


halos. This will be useful for any future numerical or semi- 
analytic project that wishes to examine the effect of flybys 
using realistic distributions. We summarize our main results 
below: 

1. Flybys occur with higher Vi„fau velocities compared to 
mergers. The mean infall velocity for flybys and merg¬ 
ers are ~ 1.8 and 1.3 the virial velocity of the primary 
halo, respectively. See Fig. [^for the probability density 
of the infall velocities. 

2. Flybys typically penetrate to about the half mass ra¬ 
dius of the primary halo, but can go as deep as < 0.01 
Rvii . However, the impact parameter is strongly depen¬ 
dent on the mass ratio; more minor flybys usually have 
smaller impact parameter. See Fig.j^for the probability 
density of flyby impact parameters. 

3. There are two classes of flybys: slow interactions that 
skirt the outer halo, and fast flybys that delve into the 
core. 

4. Flybys typically last roughly 4-5 tcmss for z < 5. There 
is a slight mass dependence - less massive primary ha¬ 
los under-go shorter-lived flybys. At high z, most flybys 
are quite short in units of crossing time. See Fig.|^for 
the distribution of flyby durations. 

5. Flybys and mergers have a similar dependence on mass- 

ratio, Major flybys are very rare (even more so 

than major mergers); the largest q we found was ~ 0.4 
in our small volume. However, at high z, flybys oc¬ 
cur between similar mass halos. At lower z, massive 
primary halos undergo fly_^s with small halos, with a 
typical q ~ 10“^. See Fig.|^ 

6. Deep encounters strip about 50% of the secondary halo 
mass, while grazing encounters exhibit negligible mass 
loss (and may even cause a gain in secondary mass). 
See Fig.|^ 

7. We find that flybys are usually highly perturbative for 

secondary halos (of order unity perturbation) but are 
typically extremely weak for the primary halos (~ lO”'^, 
see Fig. [^. Since the mass-ratio for a typical flyby, is 
^ 10“^, the perturbation, AL /E cx q^, is tiny. However, 
since massive primary halos undergo a large number of 
flyb^, some of these flybys are highly perturbative. In 
Fig. 1^ we show that highly perturbative flybys occur at 
> 1 per Gyr for all halos > for z < 4j^ 


5. CONCLUSIONS This work was conducted in part using the resources of the 

In this paper, we characterized halo flybys by their infall ve- Advanced Computing Center for Research and Education at 

locity, impact parameter, duration, mass-ratio distribution, as Vanderbilt University, Nashville, TN. We also acknowledge 

well as the mass loss and perturbation due to the flyby on both support from the NSF Career award AST-0847696. 

APPENDIX 

CALCULATING THE PERTURBATIONS EROM FLYBYS 

Flybys are high-speed encounters between halos - therefore, one could assume that the impulse approximation is valid. Flybys 
will impart an impulse on the secondary halos, resulting in an initial increase in the kinetic energy of the secondary. If the 
secondary was in equilibrium prior to the flyby, this impulse will throw it out of equilibrium and cause it to expand. Under the 

* So, contrary to standard lore, the answer to the question in the paper-title 


is not always ‘no’ 
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impulse approximation, the change in energy of the secondary (??) is; 


A£2 



/ Ml \ ^ r^ep^ 
VVsep/ 


(Al) 


For an extended perturber, the right hand side of Eqn. |^is modified with an additional factor (?), f{b), where f(b) changes 
from 1 for = oo to 0 for b ~ 0. In that modified form, the largest A £’2 occurs for rsep ^ R\. The ‘grazing’ flybys that we 
study approximately have rjgp ~ should thus have the strongest perturbations. However, both the impulse and tidal 

approximation disregard resonant effects in the secondary. To get a more accurate estimate of the overall change in the potential 
due to a flyby, we use the empirical results in Table 3 from ?. Briefly, the change in energy of a perturbed system is given by: 
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(A2) 


where, fl = ^ y 7 ^^, Rsep and Vsep refer to the relative separation and velocity of the centers of the two halos, M is the total 

mass of a halo, G is the universal gravitational constant, Ryir is the virial radius of a halo, <!> is the total potential energy of the 
halo and the subscripts 1 & 2 refer to the primary (more massive) and the secondary halo respectively. For a given pair of halos, 
AE(E is larger for the secondary halo. For penetrating encounters, we calculated the perturbation, AE /E, of the primary using 
the fit including the effects of damped modes in Table 2 of ?; 


AEi 


( —^ xKx Vsep “ where, 
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fd.0,1.96) 
1(0.5,1.90) 


0<l7'< 1, 
l<b<2. 


(A3) 


Here, b' = Rsep/Rhaif, is the impact parameter scaled in terms of the half-mass radius of the primary, and Vsep is the relative 
velocity in units of 200km/s. However, the internal and external fits are not continuous functions with impact parameter. We 
modified the fits to allow for a smooth transition between the internal fit (for b <2x primaryRhaif) and the external fit (for 
b > 2primary Rhaif): 
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exp(-0.5/3) exp(/3l7)V; 


sep 


(A4) 


whe re, w e have replaced K with exp(-0.5/3) x exp(/3l7). We used a= 1.95, /3 = -1.0. This provides a continuous fit for AE/E in 
Eqn. A4 while (roughly) preserving the parameter values from the discrete fit. 












